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ABSTRACT 

We present the results of a variety of simulations concerning the evolution of 
multiphase (inhomogeneous) hot interstellar medium (ISM) in elliptical galaxies. We 
assume the gases ejected from stars do not mix globally with the circumferential gas. 
The ejected gas components evolve separately according to their birth time, position, 
and origin. We consider cases where supernova remnants (SNRs) mix with local ISM. 
The components with high metal abundance and/or high density cool and drop out 
of the hot ISM gas faster than the other components because of their high metal 
abundance and/or density. This makes the average metal abundance of the hot ISM 
low. Furthermore, since the metal abundance of mass-loss gas decreases with radius, 
gas inflow from outer region makes the average metal abundance of the hot ISM 
smaller than that of mass-loss gas in the inner region. As gas ejection rate of stellar 
system decreases, mass fraction of mass-loss gas ejected at outer region increases 
in a galaxy. If the mixing of SNRs is ineffective, our model predicts that observed 
[Si/Fe] and [Mg/Fe] should decrease towards the galactic center because of strong iron 
emission by SNRs. In the outer region, where the cooling of time of the ISM is long, 
the selective cooling is ineffective and most of gas components remain hot. Thus, the 
metal abundance of the ISM in this region directly reflects that of the gas ejected from 
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stars. Our model shows that supernovae are not effective heating sources in the inner 
region of elliptical galaxies, because most of the energy released by them radiates. 
Therefore, cooling flow is established even if the supernova rate is high. Mixing of 
SNRs with ambient ISM makes the energy transfer between supernova explosion and 
ambient ISM more effective. 


Subject headings: galaxies : intergalactic medium - galaxies : interstellar matter - 
galaxies : X-rays 


1. INTRODUCTION 


Observations by X-ray satellites have shown that the X-ray properties of bright elliptical 
galaxies can be explained by thermal emission of hot interstellar medium (ISM). The estimated 
temperatures are around 1 keV and X-ray luminosities are typically around 10 40 erg. Theoretical 
arguments indicate that the ISM is inhomogeneous; Mathews (1990) estimated that the IMq 
of metal ejected by each supernova event into the ISM is trapped locally within the hot bubbles. 
Since in elliptical galaxies, there is no overlapping of expanding supernova remnants after galactic 
wind period (Fujita, Fukumoto, & Okoshi 1996 ; hereafter Paper I), it is expected that this 
inhomogeneity persists for a long time. The recent observation of ROSAT supports this idea 
by showing existence of components with several temperatures in elliptical galaxies. (Kim & 
Fabbiano 1995) 

Based on these arguments, we studied the evolution of the multiphase (inhomogeneous) 
ISM in Paper I. We showed that the gas components with high metal abundance and high 
density cool and drop out of the hot ISM faster than other components. As gas ejection rate 
of stellar system decreases, gas components with low metal abundance dominate in a galaxy. 

As a result, the average metal abundance which we predict is lower than that predicted by the 
previous homogeneous ISM models (e.g. Arimoto & Yoshii 1987 ; Matteucci & Tornambe 1987 
; Loewenstein Sz Mathews 1991 ; Renzini et al. 1993 ; Mihara & Takahara 1994 ; Matteucci 
Sz Gibson 1995 ; Fukumoto Sz Ikeuchi 1996). The low metal abundance is consistent with the 
observations. (Awaki et al. 1994 ; Loewenstein et al. 1994 ; Matsushita et al. 1994 ; Mushotzky 
et al. 1994 ; Kim Sz Fabbiano 1995 ; Davis Sz White 1996 ; Matsumoto et al. 1996). We also 
predicted in Paper I that the heating by supernovae is not effective because the gas components 
originated from supernovae remnants have high metal abundance, and their energy radiates before 
transfered to the circumferential ISM. 

However, since our model was a simple one-zone model, it could not predict spatial variation 
in the ISM, which is recently observed by ASCA (Mushotzky et al. 1994) and ROSAT (Kim Sz 
Fabbiano 1995 ; Rangarajan et al. 1995 ; Irwin Sz Sarazin 1996). In this paper, we investigate 
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the radial evolution of the ISM assuming that the gases ejected from stars do not mix with the 
ambient ISM globally. 

The paper is organized as follows. In §2 we describe our method for simulating the evolution 
of the multiphase hot ISM in elliptical galaxies. Then, in §3, we solve these equations numerically 
for some typical models. Our results are summarized in §4. 


2. ASSUMPTION AND BASIC EQUATIONS 


We consider the evolution of the hot ISM after galactic wind period. We assume for simplicity 
that an elliptical galaxy is spherically symmetric. The model galaxies we used here are the same 
as those of Loewenstein & Mathews (1987). The stellar density distribution is 

P*(R) = p 0 *[l + CRARa*) 2 r 3/2 (l) 


for R < Rt, where R is the distance form the center of the galaxy and Rt is the tidal radius. The 
density distribution of dark halo is 

Ph{R) = P0h\\ + {R/Rahj 2 }- 1 (2) 


for R < Rf. The subscripts a and 0 refer to core and central properties. Both densities are 
assumed to vanish for R > Rt. The properties of the model galaxies are shown in Table 1. The 
stellar velocity dispersion profile, <t*(I 2), is derived by equation (23) of Sarazin & White (1987), 
which is valid for an isotropic velocity dispersion. 


The ISM of the galaxy is divided into finite number zones, 0 < Ri(t) < ... < Rj(t) < ... < 
R m (t ), Each zone moves with the ISM. The zone, Rj-i < R < Rj, is called the j th zone. When 
Rj—i > Rt, the jth zone is deleted. Furthermore, when Rj < 0.12 kpc, the jth zone is also deleted 
to avoid too many time-steps. Each zone is assumed to be uniform on average, although it contains 
many gas components (hereafter called “phases”). For simplicity, we assume that the phases 
do not globally mix with ambient gas and comove with the zones which they belong to. These 
assumptions are valid as long as magnetic field suppresses relative motion of the phases (Hattori, 
Yoshida, & Habe 1995). We will discuss the validity in § [2.1.3 . The phases are classified by their 
origin, birth time, and zone which they belong to. We consider the three types of the origin, 
that is, gas ejected through stellar wind (mass-loss gas), shell and cavity of supernova remnants 
(SNRs). The phases originated from mass-loss gas are called “mass-loss phases” indicated by the 
index ML. Since the cavity of a SNR is buoyant, it floats and breaks part of the shell, and mixes 
with the ambient ISM locally. These partially broken shell and mixed cavity evolve as separate 
phases. They are called “bored shell phases (BSPs)” and “mixed cavity phases (MCPs)” indicated 
by the indices bs and me, respectively. We divide time into a finite number of steps for the jth 
zone, 0 < toj < t±j < ... < Lj < ... < t n j (see equation ([5(i|)), where to,j(= to) and t n j(= tf ) are 
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the start and end time of the calculation, respectively. The both times are independent of the 
zones. We define t = to(= 0.5Gyr) as the time when the galactic wind stops. In each time-step, one 
phase is born for each type. The phase of type a ( ML , bs , or me) which is born in tj-ij < t < t % j 
and in the jth zone is called the (i,j,a )-th phase and has its own temperature T^ ,a \ density 
metal abundance Z^ ,a \ and mass M^ ,a \ Each phase radiates its thermal energy and 
evolves. We assume that energy transfer between the phases is worked by pressure. Therefore, 
the phases are not independent each other. We ignore thermal conduction for simplicity. The 
phases are assumed to be in pressure equilibrium because sound crossing time of an elliptical 
galaxy is far shorter than its age. Exceptionally, the phases whose temperatures become below 
T cr it(= 10 5 K) are not considered to be in the pressure equilibrium because the condition of the 
pressure equilibrium will broken down for the phases due to the high cooling rate. They are 
assumed to cool immediately and drop out of the hot ISM. The gas left after galactic wind in the 
jth zone is also considered as a phase called “zero-phase”or (0, j, 0)-th phase. 


In § |2.1| , we describe the formation of the phases in the jth zone at t = Uj(l 
unless otherwise mentioned. We will often omit the subscripts i and j in § |2.1| . 
the evolution of the phases and the galaxy. 


< i < n, 1 < j < m) 


In 


2.2 


we describe 


2.1. FORMATION OF EACH PHASE 

2.1.1. Mass-Loss Phases 


We assume that the temperature of the gas ejected by the stellar wind immediately becomes 
equal to the stellar temperature of the galaxy T*(R) = [cr*(R)] 2 [imn /&b, where fi is the mean 
molecular weight (= 0.6), mu is the mass of the hydrogen atom, and is Boltzmann constant. 


Under this assumption, the initial temperature and density of the mass-loss phase are then 
given by 


= T^Rj) , 

Ji.i.ML) u \ hmnP b) (tij) 

P kaURj) 


( 3 ) 

( 4 ) 


where P^\t) is the average pressure of the ISM in the jth zone. The pressure P^\t) is obtained 


by solving the evolution equations of the galaxy described in §2.2. 


The gas and iron mass of the mass-loss phase at its birth time, 
M^ ML) (t hJ ), respectively, are given by 


and 


M^ j ’ ML \t itj ) = / V L 


U,ML) 


(' t)dt , 


( 5 ) 
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L 


fti-i, 


a, ml) 

Fe 


(t)dt , 


( 6 ) 


where L±' AIL ^ and L^ ML ' > are the gas and ii 
respectively. Since the time-scale of star formation in an elliptical galaxy, which is typically 10' -8 
yr, is short enough compared with the galaxy age (e.g. Arimoto & Yoshii 1987), we assume that 
the stellar system of the galaxy formed at t = 0 simultaneously. Thus, L+’ ML ' > and Lp’ ML ' 1 are 
given by 


j Fc aic unc gas and iron mass ejection rates from stars in the jth zone, 


L 


( 3 , ML) 


(■ t) = f 


dm(r ) 


dr 


<j){rn{r = t))M*\t) 


T=t 


= Z ML f 


dm(T ) 


dr 


4>{m{r = t))Mi j \t) , 


( 7 ) 


( 8 ) 


T — t 


where <f)(m) and Zml are the initial mass function (IMF) and the iron abundance of the mass-loss 
gas, respectively. The power of the IMF is taken to be 1.35. Total mass of the stars in the jth 
zone, is given by 


rRj{t) 


M ( J\t)= I 4t xR 2 Pi< {R)dR 

The relation between the mass of a star m(r), and its lifetime r, is given by 

log 10 m(r) = 1.983 — 1.054^/log 10 r + 2.52 , 


( 9 ) 


( 10 ) 


where the mass is in units of solar mass ( Mq ) and the lifetime is in units of Gyr. (Larson 1974). 
Because we consider the evolution of the hot ISM after galactic wind stops (t > to = 0.5Gyr), 
this equation implies that we do not have to consider stars with mass larger than 2.85 Mq in our 
calculation period. We assume that stars with mass in the range of 0.1 — 2.85 Mq lose their masses 
by stellar winds ; we simply assume that the mass loss occurs instantaneously at the end of the 
life given by the fraction / of the initial mass of a star with m, 


f 0 for m < 0.7 , 

/ = l 0.42m for 0.7 <m< 1.0 , (11) 

[ 0.8 — 0.43/?b for 1.0 < m < 2.85 , 

(Koppen & Arimoto 1991). 

Although the metal abundance distribution of the stars or mass-loss gas is not well-known 
beyond the effective radius, it is assumed to be 

Z ML (R) = Z ML (Re)(R/Re)- 1/2 (12) 


(Arimoto et al. 1997), where R e is the effective radius and we assume that R e = 12 R a ±. 

The density p^’ ML \tij), the temperature the mass M^ l d’ ML \tij), and the 

iron mass j) determined by above equations give the initial conditions of the evolution 

equations of the phases (see § |2.2| ) . 
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2.1.2. Bored Shell and Mixed Cavity Phases 


Next, we describe the formation of the shell and cavity phases. We consider only Type la 
supernovae (SN la) because we are concerned with the evolution of the hot ISM after galactic 
wind stops. 


In Paper I, we assumed that the SNRs consist of two parts, outer shell and inner cavity 
region, and that they evolve as separate phases. We called the former “shell phase” and the latter 
“cavity phase”, denoted by the indices s and c, respectively. When the distance from the center of 
the SNR is r and the SNR radius is r s , which is referred as “shock front radius” in Paper I, the 
inner cavity region and the outer shell region correspond to r < (1 — k)r s and (1 — k)r s < r < r s , 
respectively, where k is representing the width of the shell. The SNR radius is given by 


r s = 




11/3 


25(7 + 1)J 
pU) 


70 


10 5 cnr~ 3 K 


(P(- ? ))- 1 /3(2.02£ , sn) 1/3 

(13) 

Y 1/3 f £sn \ 1/3 

) \ 10 51 erg/ pc ’ 

(14) 


where j(= 5/3) is adiabatic constant, and Esn(= 10 51 erg) is the energy released from a supernova 
(Paper I). 


The density, temperature, and mass of the shell and cavity region before the cavity floats by 
buoyant force are 

Ps = 


rr s 

AttP sd(r)r 2 dr / / Anr 2 dr , 

k)r s J(l-k)r s 

(15) 

3 m c 

Pc 47 t( 1 — fc) 3 ?’f 

(16) 

_ 2pm n U s f{j) 

(17) 

_ 2fim H U c f{j) 

C “3 *b +1 ’ 

(18) 

rs 

m s = p s Anr dr , 

J (1 —k)r s 

(19) 

rn c = m S NR - rn s , 

(20) 


respectively, where psd(r) is the density distribution of the Sedov solution, U is the specific 
thermal energy, T^) is the average temperature of the ISM in the jth zone, and msNR is the mass 
inside the SNR radius. Their derivations are described in Paper I. Note that there is no physical 
meaning in using Sedov solution ; it is only an approximation of Fig.2 - 4 of Mathews (1990) as 
we discussed in Paper I. 
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In Paper I, we ignored floating of the hot cavity region by buoyant force and subsequent 
mixing by Rayleigh-Taylor instability (buoyant mixing) for simplicity, although Mathews (1990) 
indicated that it is effective. In this paper, we consider the effect as the followings. 

The time-scale for the cavity region to move by one diameter and to mix into the local ISM 
is <; 10 6-7 yr (Mathews 1990). Since it is less than the time-scale of the evolution of elliptical 
galaxies, we assume that the cavity region moves by one diameter and mixes into the local ISM 
immediately after its birth. Figure 1 gives an outline. The region, A + B + C, that is, the local 
ISM, part of the shell region, and the cavity region, is assumed to mix uniformly. Hereafter, we 
refer to the region A + B + C as the mixed cavity region, and refer to the region D as the bored 
shell region, denoted by the indices me and bs, respectively. They evolve as separate phases. We 
call the former “mixed cavity phases (MCPs)” and the latter “bored shell phases (BSPs)”. 

The density, temperature, and mass of the bored shell region and mixed cavity region are 


Pbs — Ps j 


( 21 ) 


Pmc — 


m c + p s V B + ptDyA 


T — - 

1 mr* — 


Pa + V B + Pc 

Tbs = T s , 

2 nmn -Esn - m bs U s 


+ T u) 


3 hn m mc 

mbs = PsV \B , 

m mc = m c + p s V b + p u) Pa , 


( 22 ) 

(23) 

(24) 

(25) 

(26) 


where Pa, Pb, and Pc are the volumes of the region A, B, and C, respectively, and p-P is the 
average density of the j’th zone. 


In Paper I, we also assumed that the metal ejected from a Type la supernova is uniformly 
dispersed within the SNR radius. However, Jun, Jones, & Norman (1996) showed that most of 
the metal is trapped in the cavity region even in inhomogeneous medium. Thus, in this paper, we 
assume that the metal ejected from a Type la supernova is trapped in the mixed cavity region. 


After the buoyant mixing, the iron mass in the bored shell and mixed cavity region are 

m Fe ,bs = m bs ZW , (27) 

TflF'e,mc — m Fe + m mc Z , (28) 

where m Fe is the iron mass ejected by a supernova and Z^T is the average iron abundance of the 
ISM occupied by the SNR. 


The densities and temperatures determined so far are adopted as the initial values of the BSP 
and MCP for the evolution equations of the phases described in § |2.2| as follows : 

p^Kuj) = Pbs , 


(29) 




(30) 

(31) 

(32) 
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T^’ bs \t id ) = T bs , 

= Pmc , 

rp(i,j,mc) (± \ _ rp 

- L \ l ^i3 ) — J -mc * 


They are determined uniquely if EgN, m pro , n^Fe, k, p^\ T^\= pm^P^ / r pC?)), and Z^ are 
given, where m pro is the mass of the progenitor star and is included by mgNR We fix Eg^, 
m pro , and m-pe by giving typical values, that is, m pro = 1.4Mq, rriF e = 0.5Mq, and ^sn = 10 51 
erg. We also £x k = 0.3. The other parameters, p^\ T^\ and Z ^ are time-dependent and 
determined by solving the evolution equations of the galaxy described in §2.2. For examples, 
for = 1.67 x 10~ 26 g cm -3 , T^) = 10 7 K, and = I.OZq, we can estimate each quantity 
: p bs = 2.4 x 10~ 26 g cnh 3 , p mc = 8.2 x 10 -27 g cnr 3 , T bs = 1.5 x 10 7 K, T mc = 1.5 x 10' K, 
m Fe,bs/ m bs = IZq, and mFe,mc/™ m c = 2.9Zq. As one can easily calculate from these values, 
the pressures of the MCP, BSP, and circumference ISM are not identical, that is, they have not 
been in pressure equilibrium yet. Therefore, we once regard equations (^) - (^) as the values 
at t = ti-i j and as the initial conditions of the iteration to derive the values at t = t l{1 (see the 
remark after equation (|56|) ). Thus, after the iteration, the MCP, BSP, and circumference ISM are 
in pressure equilibrium. Although the cooling time of the MCP is shorter than that of the cavity 
phase in Paper I, it is longer than that of the shell phase in Paper I. This means that more iron 
ejected by supernova explosion remains in hot ISM until t = tf than in Paper I. The emission from 
the MCP have influence on X-ray spectrum (see §j3i). 


The gas and iron mass of the BSP and MCP at their birth time are given by 




(33) 


[ U ' j Li j ’ mc \t)dt, 

(34) 

M^ bs \tij) = 

f Uj Lt bs) m, 

(35) 


r 4 r'w, 

^ t-i— l,j 

(36) 


where L*’ bs \ L+’ mc \ Lp^ bs \ and Tp e ’ m ^ are the mass and iron production rates of the BSP and 
MCP in the jth zone, respectively ; they are given by 


L[ J,bs \t) = r SN (t)m bs (t)Mi l] (t) , 

(37) 

L{ j,mc \t) = r S N(t)m mc (t)Mi j \t ) , 

(38) 

Lp^ bs \t) = rsN(t)m Fe , b s(t)Ml j \t) , 

(39) 

L^ e mc \t) = rsN{t)mFe,m C {t)Mi J \t) , 

(40) 


where rgrc(i) is the SN la rate per unit stellar mass. The time dependence of the SN la rate is 
assumed to be rg^(t) oc 5 (David, Forman, &; Jones 1990). The normalization is given by the 
SN la rate at t = tf, which is specified later. 
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2.1.3. Can Magnetic Field Suppress Mixing of the Phases? 


Nulsen (1986) pointed out that perturbations can be pinned in cooling flows when magnetic 
stresses are able to suppress their relative motions. This occurs when the Alfven speed, v\ is larger 
than the terminal velocity, vt defined by the force balance between ram pressure and the excess 
gravity. Hattori et al. (1995) confirm this numerically; they show that magnetic field can suppress 
relative motion in a non-linear perturbation and can help the growth of thermal instability as long 
as the size of perturbation is smaller than the critical value 


A 


crit 





(41) 


where (3 is the ratio of the gas pressure to the magnetic pressure, and L is the scale height of 
pressure. The density contrast, 5, is defined as 


§= P^Pf 

Pb 


(42) 


where p and pb are the density of the perturbation and unperturbed flow, respectively. 

For the case of a temperature of 10' K and a density of pb = 1.67 x 10 -26 gem -3 , and there 
is a 1 pG magnetic field, one finds that (3 = 350 and r s = 70 pc ; in the following argument in 
this subsection, we will use these parameters and take L = 5 kpc. When k = 0.3, spatial scale 
lengths of cavity and shell region respectively are 50 pc and 20/2 = 10 pc, following the definition 
of Hattori et al. (1995). Since the initial value of 6 for the shell region is less than one, the critical 
wavelength is A cr it > 15 pc. This indicates that the shell region can be supported by magnetic 
field before the cooling becomes efficient. On the contrary, one can show that the cavity region, 
before mixing, cannot be supported, if equation (|a]) can be applied to less dense perturbation for 
151 ~ 1 ; it will move upwards and mix with the ambient medium as we assumed in §|2.1.2|. 


As each phase cools, the density contrast increases. Hattori & Habe (1990) show that if 
magnetic tension supports a perturbation until 5 > 10, the gas in the perturbation can cool to 
10 4 K. From equation (41), one finds A cr it = 1.5 pc when 5 = 10. Since the spatial scale length 
of gas blobs composing a mass-loss phase is typically ~ 1 pc (Mathews 1990), the blobs can cool 
to 10 4 K without mixing and drop out of hot ISM. On the other hand, since the spatial scale 
lengths of mixed cavity region and bored shell region are 50 and 10 pc, respectively, they cannot 
cool to 10 4 K without mixing unless they fragment into small pieces before the cooling becomes 
effective (^ 10 7-8 yr). However, Jun et al. (1996) show that a SNR in inhomogeneous medium 
has complicated structure ; thus, we do not think that each region evolves as one gas blob for 
a long time. Moreover, the estimation of Hattori & Habe (1990) does not include the effect of 
metal abundance fluctuation which makes gas blobs with higher metal abundance cool faster, and 
makes the value of 5 until which magnetic tension should support them smaller. Furthermore, 
since a SNR amplifies magnetic field (Jun &; Norman 1996), A cr ;t may be large around it. Thus, 
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we assume in §J3| that the MCP and BSP also comove with ambient ISM, although we consider an 
additional model in which the MCP and BSP mix with ambient ISM in wider region. 

The above estimations may be too simple ; strength variation and structure of magnetic 
fields may affect the evolutions of gas blobs. However, since a full treatment of convection in 
inhomogeneous medium would be very difficult, we consider comoving flows as a first-step of the 
research. 


2.2. EVOLUTION OF THE PHASES AND THE GALAXY 


The energy equation for the th phase for t > Uj is given by 

d (k B T^’ a \t)\ k B T^ a \t) d 


7 — 1 dt 

,(m» 


dm, H 


prn B dt 


p^ a \t) = ~(n Rj ’ a) ) 2 A{Z Rj ’ a \T (i ' j ’ a \t)) , (43) 


where rv e ; is the electron density of the (i,j,a )~th phase and A is the cooling function 
approximated by 


2.1 x 10“ 27 1 + 0.1 


\ / y(U» ' 


0.5 


Z© 


K 


+ 8.0 x 10 -17 0.04 + 


(ergs cm 3 s 7 ) 




-l.o 


Z© 


K 


(44) 


This is an empirical formula derived by fitting to the cooling curves calculated by Bohringer & 
Hensler (1989). 


In our model, we assume that all existing phases in a certain zone have same pressure. Thus, 


for tij < tij, 


d m n 


(45) 


We assume that the ISM is in hydrostatic equilibrium, that is, 

ApU) _ d) GM ( R j) 

ARj ~ P ' ( Rj ) 2 ’ 


(46) 




p 4 j \t u ) 


(47) 


where Mg' 1 and are respectively the total gas mass and volume of the jth zone, 

A = pC+i) _ p(J) f and A Rj = Rj — Rj-i- We define as the pressure just outside Rj-i- 
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The gravitational mass within R is given by 

M(R)= [ R 47rR 2 (p* + p h )dR. 
Jo 


(48) 


The total gas mass, the iron mass, and the gas volume of the j the zone are the summation 
of those of each phase ; 

M g (i) W = E 51 (49) 

a “)>T crit 

M SeW = E E (50) 

“ i, T(''J'“)>T cri t 

^■ ) w = E E (si) 

“ i, T( i 'J'“)>T crit 

where 

K w '“) (i) = (t)/ p 0d,«) (*) . (52) 

As mentioned above, the phases whose temperatures are below T cr ; t (= 10' 5 K) are not included in 
the summation. 


Note that part of the SNRs is composed of pre-existing phases ; the masses of the pre-existing 
phases are reduced by the occupation by the SNRs. Thus, the mass and iron mass of the 
pre-existing (i,j,a )-th phase at t = tij (tij < tij ) are given by 




t=tl-1 , 


«£**’(««) = - Kmj 


(h3,a) i 


Aj) 


V&(t) 


(53) 


(54) 


t— 


respectively, where Mgjj RZ is the mass occupied by SNRs during ij <t< tij, and is given by 


Af« 


(j) _ 


SNR,/ 


^l — l,j 


(t)V s rsa(t)M*\t)dt , 


(55) 


where Vs = Va + Vb + Vc + VR- 


We take the sound crossing time as the time-steps for the calculations, 


ti+i,j ti,j 


ARj(tij) 


7 


ksfW 


- 1/2 


( 56 ) 


We derive pb>J- Q )(/ /j -) and T^ J - a Hti j) from p(*d>“)(^_ Rj .) an d TV J > Q )(/;_] .,■) by iterating 
Eqs.(^3|), (^) , (p46l) , (^9|) , (fill) , and (|53|) until they converge. 







3. MODEL RESULTS 


In this section, we present results of the evolution equations described above. We reduce 
the number of the free parameters by giving typical fixed values to some of them. First, for the 
parameters regarding the whole galaxy, we take to = 0-5 Gyr, tf = 10 Gyr. 

We assume that the gas distribution and abundance distribution at t = to are 

Pgas{R ) = POgasY + {R/Ra*Y\ > (57) 

Z 0 (R) = Zu{Re){R/R e )- 112 (58) 

for R < R t , respectively. The normalization, po gas , is determined by assuming that the total gas 
mass at t = to is 1/30 of the total stellar mass. The normalization, Zo(R e ), is a parameter. We 
give the initial zone radius so that enough spatial resolution is obtained at t = tf. The initial 
number of the zone, m, is 30. We can write 

= V&(to) = y [Rj(to) 3 - Rj-l(to) 3 ] , (59) 

pRj (to) 

M(°^°)(fo) = M^(to) = / 4irR 2 p gas (R)dR , (60) 

8 jRj-Rto) 

M^’ 0) (to) = Mg g (t 0 ). (61) 

We take the temperatures of the zero-phases 

T^°\to) = T i ,(R j (t 0 )). (62) 

We solve the basic equations described in §j| for the galaxy models whose details are given 
in Table 1. The SN la rate is normalized by the values at t = tf(= lOGyr) shown in Table 1. 
The values are expressed in units of SNu, that is, the number of SNe per 10 10 /U 2 L b q per 100 yr 
(Hq = 100/ikm s -1 Mpc -1 ; we set h = 0.5). The ratio of the mass to luminosity is taken to be 
8 Mq/Lq, and the rate for SN la is normalized by using this ratio. Models A3-A6 and B3-B6 
are calculated to see the effect of metal-abundance distribution for the mass-loss gas. In these 
models, the mass-loss phase born in each time-step is divided into two phases with equal masses 
and different abundances corresponding to the two figures for ZMh(Re) (Table 1). On the other 
hand, in models A1 and A2, each component has the same abundance or the mass-loss gas has one 
kind of the abundance. In model A5, shell and cavity regions mix into circumference ISM of 26 
times their volumes, respectively. This means that shell and cavity region mix into circumference 
ISM in wider region than other models. 

Figure 2(a) shows the time evolution of Rj (1 < j < m) in models A1 and A2 ; this figure 
shows the influence of SN la on the gas evolution. Cooling flow is established in the inner region. 
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Note that in other models, the cooling flow is also established. Although the supernova rate of 
model A2 is larger than that of model Al, the inward velocity in model A2 is not reduced. This 
fact indicates that supernovae are not effective heating sources of the ISM in the inner region 
of the galaxy. The reason is that the BCP and MCP cool effectively owing to their high metal 
abundance (see §2.1.2), and that they have large thermal energies when they are born ; the sum 
of their energy is 

3 kf 


Ebs + E mc — PlgN + 


2/rmn 


pVs > -Esn 


(63) 


except for models A5 and B5. This indicates that more energy than that supplied by supernova 
explosion radiates from the BSP and MCP in a short time. Thus, supernovae cannot heat the 
ISM effectively. In outer region, the outflow velocity in model A2 is larger that in model Al ; 
supernovae are effective as heating sources. In this region, cooling times of the BSP and MCP 
are so long that they do not radiate their thermal energy so much. Thus, the energy released by 
supernovae is transfered to the circumferential ISM before the cooling. 


In Figure 2(b), outward velocity of model A5 is larger than that of model A3. This is because 
mixing of SNRs makes their cooling time longer. Therefore, the energy transition between SNRs 
and circumferential ISM is more effective in model A5. 


Figures 3-5 show the distributions of the density, luminosity averaged temperature, and 
metal abundance at t = tf (at t = 9.4 Gyr for model A5 because the radius of the innermost 
zone is ~ 20 kpc at t = tf). We confirmed that, for models in Paper I, the luminosity averaged 
temperatures are almost same as those derived by spectral simulations which we did in Paper I. 
On the contrary, the luminosity averaged metal abundances are not consistent with those derived 
by the spectral simulations. Thus in this section, we simply use Z m = as the metal 

abundance, and later, derive X-ray spectrum of the model galaxies. The units of metal abundance 
is the solar abundance Zq, 1.7 x 1CP 3 . Since our model does not have enough spatial resolution 
in the central region of the galaxy, we infer the distributions of the density, luminosity averaged 
temperature, and metal abundance by presenting the values of the zones which disappear just 
before t = tf. The values shown in Table 2 are the ones when the zone radii are ~ 0.3 kpc. For 
model A2, we stopped calculating the evolution of the innermost zone and deleted it before its 
radius decreases to 0.3 kpc. This is because masses of some phases in the innermost zone become 
negative by equation (|53|), because of large Mgj] R l and density contrast among phases. In order to 
overcome this defect of our model, the exact treatment of the evolution of supernova remnants in 
inhomogeneous medium is required. However, we think this deletion does not affect the evolution 
of the outer zones, because when we stop the calculation, the volume of the innermost zone is far 
smaller than those of the outer zones. 

Figure 3 show that the density varies crudely as p oc R~ L5 for R ^ 1 kpc. The distributions 
of density are almost independent of the models. This reflects that the density distributions are 
mainly decided by energy injection rate and vary as oc pj (Sarazin, & White 1987). The high 
central density of models A5 and B5 reflects small density contrast among the phases in the 
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innermost zones (Table 2). 

The temperature gradually rises toward the galactic center (Figure 4). This temperature rise 
is caused by adiabatic heating owing to the gas inflow (Figure 2) ; in our models, the phases with 
low-cooling rates are brought to the center and compressed, because high-cooling phases drop 
out and the volume of the inner zones reduces. In the central region, the temperature falls down 
(Table 2). In this region, all phases start to cool because of high pressure. 


The metal abundances, Z m , are very different among the models A1-A3 (Figure 5(a)). We 
show the metal abundance distribution of mass-loss gas when Zml(^) = 1 Zq (equation (HD) in 
Figure 5 for reference. In model A2, the high supernova rate makes the metal abundance high. 
This high metal abundance is inconsistent with the observations (Awaki et al. 1994 ; Loewenstein 
et al. 1994 ; Matsushita et al. 1994 ; Mushotzky et al. 1994 ; Kim Sz Fabbiano 1995 ; Davis 
Sz White 1996 ; Matsumoto et al. 1996). This may indicate that the real supernova rate is 
lower than that derived by the observation (~ 0.2 SNu ; Tammann 1982). In the central region, 
high-abundance phases generally cool and drop out of hot ISM faster than other phases, although 
it depends on their temperatures and densities. This makes the average metal abundance in that 
region small. However, since the negative metal abundance gradient of the stars or mass-loss gas 
cancels this effect, the metal abundance of the hot ISM is constant or rises toward the center 
of the model galaxy except for model A6 and B6 (Figure 5). Since this result depends on the 
assumption of equation (12), the observation of the metal abundance of the stars beyond the 
effective radius is encouraged. Moreover, gas flows inward from the outer region where the metal 
abundance of mass-loss gas is small. This means that as gas ejection rate of the stellar system 
decreases, mass fraction of low-abundance phases born in the outer region increases in a galaxy. 
In the outer region, the selective cooling is not effective and the metal abundances of the hot ISM 
directly reflects the metal injection from stars including supernovae. 


The metal abundance in models A3 is smaller than those in models Al (Figure 5(a)), 
because the mass-loss phases with higher metal abundance cool faster than those with lower metal 
abundance. Compared with model A3, model A4 shows higher abundance (Figure 5(a) and 5(b)). 
In model A4, part of high abundance gas of the zero-phase is occupied by SNRs which survive 
until t = tf. In model A6 and B6 , mass-loss phases with low metal abundance dominate in the 
central region at t = tf. 

Although the metal abundances of the models except for A2 is low, they cannot be compared 
directly with the ASCA results (^ 0.5Zq). Since the density of phases with high abundance are 
generally high, emission from them can affect X-ray spectrum. Therefore, we simulate the spectra 
at t = tf, as we did in Paper I, for models A3 - A6. We simulate the spectra of innermost zones 
and those of zones at ~ 20 kpc. Although we do not consider projection effect, contamination of 
outer zones is less than 20% for the innermost zones. For ^ 20 kpc, the variation among zones are 
little except for densities. The distance to the model galaxies is assumed to be 15 Mpc. We use 
XSPEC package (version 9.00) and response function of ASCA Solid-State Imaging Spectrometer 
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(SIS). The model spectrum are fitted by Meka model (Mewe, Gronenschild, and van den Oord 
1985 ; Mewe, Lemen, and van den Oord 1986 ; Kaastra 1992). The method is described in Paper I. 

The results are shown in Figure 6. Because of the low luminosity, Z sp at ~ 20 kpc has 
an error of ~ OAZ q. We present the spectrum of the innermost zone of model A4 as an 
example (Figure 7). The central metal abundance derived by spectral fitting, Z sp is larger 
than Z m for models A3 and A4. This is because that MCPs affect the X-ray spectrum by 
their strong metal emission. The number of MCP which can survive until t = tf is small 
compared with those of other phases; for example, for the innermost zone in model A4 the ratio 
is MLL : MLH : BSP : MCP = 161 : 14 : 146 : 97, where MLL and MLH refer to mass-loss 
phase with low and high metal abundance, respectively. On the contrary, the luminosity ratio is 
MLL : MLH : BSP : MCP ~ 2 : 1 : 4 : 4. Note that for the zone at ~ 20 kpc, the luminosity 
ratio is MLL : MLH : BSP : MCP ~ 2 : 3 : 2 : 2. In model A5, Z sp ~ Z rn in the galactic center, 
because of small abundance fluctuation among the phases. However, the luminosity of MCP and 
BSP dominates that of mass-loss phase as models A3 and A4. 

In the central region, the derived metal abundances, Z sp , are larger than the ASCA results 
(Mushotzky et al. 1994) except for model A 6 , although Z sp Zml (Figure 6 ). This may indicate 
that our model is too simple to predict metal abundance quantitatively. Alternatively, supernova 
explosion rate in elliptical galaxies may be < 0.01 SNu ; the central abundance deficit of NGC 
4472 (Irwin & Sarazin 1996) may correspond to that in models A 6 and B 6 . As discussed above, 
our model predicts that in a central region of a galaxy, iron line emission from SNRs of Type la 
supernovae, if any, should be prominent. That is, observed [Si/Fe] and [Mg/Fe] should decrease 
towards the center. 

The X-ray luminosities of the model galaxies are 10 41 42 erg s 1 for t ~ 10 10 yr ; the 
luminosities fluctuate because of the small number of the zone. 


4. SUMMARY AND CONCLUSIONS 


In this paper, we have presented a model of the evolution of the hot gas for elliptical galaxies 
after galactic wind period under the assumption that the gases ejected from stars do not mix with 
the circumferential ISM globally. The ejected gases evolve separately according to their birth 
time, position, and origin. We considered three origins of the ejected gas, that is, shell and cavity 
of supernova remnants and mass-loss gas. Furthermore, we considered the floating of the cavity 
and subsequent mixing by Rayleigh-Taylor instability. 

The main results and conclusions can be summarized as follows: 


1. The model predicts that the supernovae are not effective as heating sources of the ISM in the 
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inner region of galaxies after the galactic wind stops. In the inner region, supernova remnant 
can cool rapidly because of their high density and/or metal abundance. Since the remnants 
initially have large thermal energy, the energy ejected by supernova explosions is radiated 
and supernovae do not heat up the ISM. Thus, cooling flow is established even if supernova 
rate is large. In the outer region of the galaxies, the cooling time of the remnants is long. 
Thus, most of energy ejected by supernova explosions is not radiated and it is transfered 
into the circumferential ISM. Mixing of SNRs with ambient ISM makes this transfer more 
effective. 

2. In a inner region of a galaxy, the present iron abundance of the hot ISM can be less than 
that of the mass-loss gas or stars if the supernova rate is small, because the phases with 
higher metal abundance generally cool faster and gas inflows from outer region where the 
metal abundance of the mass-loss gas is small. However, the spectral simulations show that 
predicted metal abundances are still larger than the ones observed by ASCA in the central 
region, if the present supernova rate is ^ 0.01 SNu. In the outer region where the selective 
cooling is ineffective, metal abundance of the ISM directly reflects that of the gas ejected 
from stars. Our model predicts that iron line emission by SNRs is prominent in the central 
region, and that [Si/Fe] and [Mg/Fe] decrease towards the galactic center. 
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Note added in proof 

A recent ASCA result shows that the metal abundances in central regions of elliptical galaxies 
are ~ 1 Zq (Matsushita 1997). 


This preprint was prepared with the AAS DT^X macros v4.0. 


Figure Captions 


Fig. 1.— The schematic figure of a supernova remnant. 

Fig. 2.— The evolutions of zones for (a) model A1 (solid line) and A2 (dotted line) (b) model A3 
(solid line) and A5 (dotted line). 

Fig. 3.- The present distribution of density, (a) A1 - A3, (b) A4 - A6, and (c) B4 - B6. 

Fig. 4.— The present distribution of temperature, (a) A1 - A3, (b) A4 - A6, and (c) B4 - B6. 

Fig. 5.— The present distribution of abundance (Z rn ). (a) A1 - A3, (b) A4 - A6, and (c) B4 - B6. 
The long-dashed line shows the metal abundance of the mass-loss gas when Zy\_h{Re) = 1-0. 

Fig. 6.— The present distribution of abundance ( Z sp ) for models A3 - A6. The long-dashed line 
shows the metal abundance of the mass-loss gas when Zml(-Rc) = 1 -0. 

Fig. 7.— Simulated X-ray spectra observed with the ASCA SIS for model A4. The line shows the 
best fitting (two Meka plasma + absorption column density). 



Table 1. 


Properties of the Model Galaxies 


Model 

R a * 

(kpc) 

Po* 

(10~ 21 g cnr 3 ) 

Rah 

(kpc) 

POh 

(10~ 23 g cnr 3 ) 

Re 

(kpc) 

Rt 

(kpc) 

SN rate 
(SNu) 

^ML {Re) 

(Zq) 

Zo(i? e ) 

(Zq) 

Al 

0.50 

7.20 

5.0 

1.75 

6.0 

100.0 

0.01 

1.0 

1.0 

A2 

0.50 

7.20 

5.0 

1.75 

6.0 

100.0 

0.2 

1.0 

1.0 

A3/A5 

0.50 

7.20 

5.0 

1.75 

6.0 

100.0 

0.01 

0.2,1.8 

1.0 

A4 

0.50 

7.20 

5.0 

1.75 

6.0 

100.0 

0.01 

0.2,1.8 

4.0 

A6 

0.50 

7.20 

5.0 

1.75 

6.0 

100.0 

0 

0.2,1.8 

1.0 

B4 

0.25 

21.8 

2.5 

5.30 

3.0 

50.0 

0.01 

0.2,1.8 

4.0 

B5 

0.25 

21.8 

2.5 

5.30 

3.0 

50.0 

0.01 

0.2,1.8 

1.0 

B6 

0.25 

21.8 

2.5 

5.30 

3.0 

50.0 

0 

0.2,1.8 

1.0 
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Table 2. Properties of the Central Region 


Model 

n 

( cm -3 ) 

kT 

(keV) 

Zm 

(%) 

Model 

n 

( cm -3 ) 

kT 

(keV) 

Z m 

(Zq) 

Al 

0.2 

0.8 

2 

B4 

0.4 

1.1 

1 

A2 




B5 

10 

1.2 

1 

A3 

0.3 

0.8 

1 

B6 

2 

0.8 

0.1 

A4 

0.4 

0.9 

1 





A5 

10 

1.0 

1 





A6 

0.7 

0.6 

0.1 





















n(cm' 3 ) 








































